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Abstract 



. The well-known TE/TM decomposition of time-harmonic electromagnetic fields 

in uniaxial anisotropic media is generalized in terms of four-dimensional differential- 
form formalism by requiring that the field two-form satisfies an orthogonality condi- 
tion with respect to two given bivectors. Conditions for the electromagnetic medium 
I in which such a decomposition is possible are derived and found to define three 

I subclasses of media. It is shown that the previously known classes of generalized 

CN ■ Q-media and generalized P-media are particular cases of the proposed decompos- 

able media (DCM) associated to a quadratic equation for the medium dyadic. As 
a novel solution, another class of special decomposable media (SDCM) is defined 
by a linear dyadic equation. The paper further discusses the properties of medium 
dyadics and plane-wave propagation in all the identified cases of DCM and SDCM. 

> : 

^ 1 Introduction 

1.1 TE/TM decomposition 

The most general linear electromagnetic medium (bi-anisotropic medium) can be ex- 
pressed in terms of four medium dyadics in the three-dimensional Gibbsian vector ("en- 
gineering") representation as [1, 2] 




(1) 



where the four field vectors are elements of the vector space Ei. The number of free 
parameters is 4 x 9 = 36 in the general case. It it well known that, in a uniaxial anisotropic 
medium defined by medium dyadics of the form 

% = et{u^u^ + UyUy) + e^u^u^, = 0, (2) 
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jig = //t(U^U^ + WyWy) + lIz^Z^Z, Cs = (3) 

and satisfying 

etfJ'z - fJ'tez 7^ 0, (4) 

any time-harmonic field with time dependence e^'^* can be uniquely decomposed in two 
parts, 

E = Ete + ^tm, H = Hte + Htm, (5) 

satisfying 

■ Ete = 0, • Htm = 0. (6) 

This property was probably first published by Clemmow in 1963 [3]. In the case etfiz — 
l^'t^z = the decomposition can still be made but it is not unique. Actually, such a medium 
can be transformed to an isotropic medium through a suitable affine transformation [2]. 
The TE/TM decomposition in isotropic media has a longer history [4]. 

TE/TM decomposition in the uniaxial medium can be simply demonstrated for a plane 
wave. In fact, because the fields of a plane wave in any medium satisfy E • B = and 
H • D = 0, they also satisfy 

etE B-HtU T>^ {etfXz - l^tez){^z ■ E)(u, • H) = 0, (7) 

in the uniaxial medium (2), (3). Thus, assuming (4), the plane wave must be either a TE 
wave or a TM wave with respect to the axial direction defined by the unit vector u^. Since 
any electromagnetic field outside the source region can be decomposed in a spectrum of 
plane waves, each component of which is cither a TE wave or a TM wave, the field can 
be uniquely decomposed in TE and TM parts. The same principle is valid also for more 
general decompositions. Thus, it is sufficient to consider the decomposition problem for 
plane waves only. 



1.2 Generalized decomposition 

The TE/TM decomposition theory has been generalized to media where the fields can be 
decomposed to two parts satisfying either a-E = Oorb-H = where a and b are two 
given vectors [5]. Even more generally, [6] analyzes the occurrence of ai ■ E + a2 ■ H = 
or bi • E + b2 • H = where ai • • • b2 are four given vectors. This last decomposition was 
shown to be possible in bi-anisotropic media with Gibbsian medium dyadics of the form 



^9 = 


;^(-B^ + a2bi+b2ai), 
It 


(8) 


I- 


X X T + ;^(a2b2 + b2a2), 
It 


(9) 


I- 


z X T + ;^(aibi + biai). 


(10) 




= ^(B + aib2 + bia2), 
Zr) 


(11) 



where x, z are arbitrary vectors, r, r] are arbitrary scalars and B is an arbitrary dyadic. 
Media defined by (8) - (11) have been called decomposable media [6]. The paper [7] 
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demonstrates that yet another scalar parameter a can be added to the definitions of 
the medium dyadics. In this case the definitions (8) - (11) are generahzed (after some 
manipulations) to the form 



% = a{z X T + aibi + biBi) + r]{-W + agbi + bgai), (12) 

= ri(x X T + a2b2 + b2a2) + a(B + aib2 + bia2), (13) 

= r(z X T + aibi + biai) - a{-W + a2bi + b2ai), (14) 

Jig = -q;(x X T + a2b2 + b2a2) + r(B + aib2 + bia2). (15) 



1.3 Four- dimensional formalism 

Remarkably, the conditions (12) - (15) of the decomposable medium can be formulated in 
a compact way by applying the four-dimensional differential-form formalism. The present 
paper thus starts directly from the four-dimensional definition of decomposable media. 

It is well known that the Maxwell equations, 

dA$ = 0, dA* = 7, (16) 

and the medium equation 

* = M|<&, (17) 

have a simpler appearance in the four-dimensional differential- form representation in com- 
parison with the three-dimensional Gibbsian vector formalism [8, 9, 10]. Here the elec- 
tromagnetic fields are represented by two-forms elements of the space F2, whose 
expansions in terms of three-dimensional two-forms B, D and one-forms E, H are 

$ = B + EAe4, * = D-HAe4. (18) 

Here, 7 e F3 is the source three-form, 

7 = ^ - J A £4, (19) 

consisting of three-dimensional charge three-form q and current two-form J. In the basis 
of one-forms £», i = 1 ■ ■ ■ 4, £4 denotes the temporal basis one-form. The medium dyadic 
M e F2IE2 maps two-forms to two-forms and corresponds to a 6 x 6 matrix. It is often 
simpler to consider the modified medium dyadic Mg G E2E2 defined by 

Mg = ejvLM, (20) 

in terms of a quadrivector ejv G E4. The modified medium dyadic maps two-forms to 
bivectors, elements of the space E2. Summary of definitions and operational rules for 
differential forms, multivectors and dyadics applied in this paper can be found in the 
appendices of [11, 12, 13] and, more extensively, in the book [10]. 
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The medium dyadic M can be expanded in four three-dimensional dyadic components by 
separating terms involving the temporal vector 64, temporal one-form £4, or both, as 

M = ^-Ff' Ae4-Fe4 A^~^-Fe4 A^Ae4, (21) 
which corresponds to the matrix representation 



h;-V5- ?;'VE '''' 

One can represent the modified medium dyadic by using Gibbsian medium dyadics as 

-(ei23 + 64 A !j|p;'|(IJei23 -lg^ 64), (23) 
which corresponds to the matrix representation 



The matrix is the same as in the expression (1) which involves Gibbsian vectors. 
The two sets of 3D medium dyadics have the relations [10] 

e ^ei2z\_{%-lgW-%), (25) 
P = £i23LPa> = Pg^Jei23, (26) 

^ = ei23ei23LL(!,P;'), 1 = -Wg%- (27) 

The four-dimensional formalism allows simple definition of important classes of electro- 
magnetic media. For example, if the modified medium dyadic can be expressed in terms 
of the double-wedge square of some dyadic Q e EiEi (which need not be symmetric) as 

Mg^M]^^M=MW^\ (28) 

the corresponding three-dimensional Gibbsian dyadics satisfy relations of the form [14, 10] 

f, + apj = 0, 1^ = -!^, l^g=lg^ (29) 

for some scalar a. Thus, Zg and /I J are multiples of the same dyadic while ^g and Cg may 
be any antisymmetric dyadics. In (29), a medium defined by (28) was called a Q- medium 

for brevity. Such a medium is known to have the property of being non-birefringent 
to propagating waves. Thus, media in this class can be conceived as generalizations 
of isotropic media. Also, the class known as transformation media [16, 17, 18] largely 
coincides with the class of Q-media with a symmetric dyadic Q. 

Generalizing the definition (29) by adding a term 

Mg = + AB, (30) 
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where A, B are two bivectors, leads to the class of generalized Q-media, defined in [19]. 
One can show that, for this kind of media, the three-dimensional medium dyadics must be 
of the form (8) - (11), i.e., any generalized Q-medium is actually a decomposable medium. 
However, this cannot be inverted, because the more general set of conditions (12) - (15) 
for q; 7^ cannot be achieved with medium dyadics of the form (30). At this stage it is 
not obvious how to generahze (30) to correspond to the conditions (12) - (15). 

1.4 Hehl-Obukhov decomposition 

In many applications a decomposition of the medium dyadic based on its symmetry prop- 
erties is often useful. As shown by Hehl and Obukhov [9] (following the symmetry proper- 
ties discussed by Post [20] ) , the most general medium dyadic can be uniquely decomposed 
in three components as [9], 

M = Ml + M2 + M3, (31) 

called principal (1), skewon (2) and axion (3) parts of M. The axion part M3 is a multiple 

of the unit dyadic T^^^-^ mapping any two-form to itself while the other two parts are 
trace free. The skewon part is defined so that the corresponding modified medium dyadic 

ejv[M2 G E2E2 is antisymmetric, while the principal part Mi is trace free and e^vL^i is 
symmetric. The decomposition (31) motivates an intuitive nomenclature; for example a 

medium defined by M = Mi is called a principal medium and one with M = M2 -|- M3 is 
called a skewon-axion medium. 

2 Decomposable medium (DCM) 
2.1 Plane- wave conditions 

Let us now formulate the decomposition property in terms of four- dimensional quantities. 
Assuming plane-wave fields 

*(x) = *e^l^, *(x) = ^e^l'', (32) 

the Maxwell equations (16) for the wave one- form v and electromagnetic amplitude two- 
forms ^ become 

i/A* = 0, i/A* = 0. (33) 
These imply the following representations in terms of potential one-forms (/>, ■0, 

^ = uA(f), *=:i/Ai/?. (34) 

Thus, the electromagnetic two-forms of any plane wave satisfy the orthogonality conditions 

*A<& = 0, #A* = 0, *A* = (35) 

in any medium. Defining the dot product between two two- forms ^' or two bivectors 
A,B as 

* • * = #|(ejvL*) =eAr|(* A*), (36) 
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A • B = A|(e^LB) = £iv|(A A B), (37) 
the four-form conditions (35) can be expressed as the scalar conditions 

$•# = 0, $-* = $-M|$ = 0, *-* = *|M^-M|$ = 0. (38) 

Thus, for any medium dyadic M, the two-form $ of any plane wave satisfies a condition 

of the form _ 

*|(Q;eArJ^^^'' + /5ejvLM + 7M^- M)|* = 0, (39) 

for arbitrary scalars a, /5, 7. For the modified medium dyadic this condition becomes 

*|(aeivLT^'^^ + + 7^ • Mg)|* = 0, (40) 

because we have 

• M = M^I(ejvLM) = {£N[Mgf\Mg 

= MJ|(£^LMp)=MJ-M,. (41) 
2.2 Condition for the medium dyadic 

Let us now assume that, given two bivcctors A.B G E2, the medium has the property 
that any plane wave satisfies either A|$ = or B|$ = 0. Such waves can be respectively 
called A- waves and B- waves and the medium, in analogy with the medium defined by (8) 
— (11), can be called by the general name decomposable medium. Thus, any plane wave 
in such a medium is required to satisfy 

(#|A)(B|*) = *|(AB)|c& = 0, (42) 

for two given bivectors A, B. Following [6] let us now define the class of decomposable 
media by combining (40) and (42) and requiring that the following condition be satisfied 
for all two-forms ^: 

*|(ae^[P)^ + ^Mg + 7MJ -MJI* = *|(AB)|*. (43) 

Since the left-hand side is zero for all media, this warrants that (42) is satisfied when the 
medium is such that (43) is satisfied. Requiring that this be valid for any two-forms $ 
implies that the symmetric parts of the dyadics in brackets on both sides must be the 
same. Redefining the coefficients we can write the condition in the form 

aejv J^'^^ + /^(M^ + MJ) + 7MJ ■ = AB + BA. (44) 

Although (44) is obviously enough to define a large class of media, it is not enough for 
claiming that this class covers all media for which the decomposition condition (42) is 
satisfied. The latter question must be left for the topic of further research. 

To find solutions for the condition (44), we must separate two cases: either 7 7^ or 
7 = 0. 

• The case 7 7^ requires solving a symmetric quadratic dyadic equation and the 
corresponding class of media can be called that of (proper) decomposable media or 
DCM for brevity 

• In the case 7 = the quadratic dyadic equation is reduced to one of the first order. 
It will require separate consideration and actually defines a distinct class of media 
called that of special decomposable media or SDCM. 
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3 Solutions to the medium conditions 
3.1 DCM 

Assuming 7 7^ in (44) we can set 7 = 1 without losing generality. In this case the DCM 
condition (44) can be expressed in compact form as 

Mf -M^^aWLl^^^^, (45) 

by defining 

M; = M, + /3e;vJ^')^-DC, (46) 

C = A, (47) 

D-Mg + /3D-^(D -0)0 = 6, (48) 

a' = - a. (49) 

The condition (45) is a quadratic dyadic equation, whose solutions are derived in the 
Appendix (cf. (122)). Accordingly, two subclasses of DCM are obtained. The first class 

assumes that there exist a dyadic Q G EiEi sTich that is a multiple of Q*^^). The second 

class assumes that there exist a dyadic P G FiEi so that is a multiple of eTvLP*-^-*. Let 
us consider these two cases separately and respectively call them QDCM and PDCM. 
This nomenclature is chosen in the light of the obvious relation with Q-media [14] and 
P-media [15]. 

Redefining a, the QDCM solution of (45) as obtained from (46) - (49) must be of the 
general form 

or _ _ 

M = q;T(2)^ + Men [Q^'^ + £iv [DC, (51) 

for some normalized dyadic Q, bivectors C,D and scalars M, ct. Thus, the definition (48) 
is more explicitly 

B = mW^^ ■T> + -ip-'D)C. (52) 

Given the bivectors A = C and B one can easily solve (52) for D. It is also easy to verify 
that plane- wave fields in a medium defined by (51) satisfy the decomposition property 
(42) (see Section IVA). 

The three-dimensional Gibbsian medium dyadics corresponding to the general QDCM 
can be expressed in the general form (12) - (15). In fact, adding an axion term with 
q; 7^ to the medium dyadic of the generalized Q-medium (30) analyzed in [19], it can be 
shown to correspond to the more general class (50) of decomposable media (12) - (15). 

The second possibihty in (45), that of PDCM, yields the following forms for the decom- 
posable medium, 

Mg = aeN [T^^^^ + Mejv [^^^^ + DC, (53) 
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or _ _ 

M = aP^^ + MP(2) + En [DC, (54) 

for some normalized dyadic P, bivectors C, D and scalars M, a. The definition (48) in 
this case is 

B = MD|F(2) + 1(D • D)C. (55) 

Setting again C = A, one can easily solve (55) for D in terms of given A and B. For a — 
the medium coincides with one called the generalized P-medium whose basic properties 
have been studied in in [15]. Thus, the PDCM solution coincides with the generalized 
P-medium extended by an arbitrary axion component. 

Based on the expansions of the dyadic P and the bivector product DC in 3D components 
as _ _ 

P = P^ + 7re4 + + ^£464, (56) 

DC = (di A d2 + da A e4)(ci A C2 + C3 A 64), (57) 

the 3D components of the medium dyadic M expressed as (21) in the PDCM take the 
form [15] 

^ = MFf + aU""^^ + (£123 Lda) (ci A C2) , (58) 

f' = -TT A + (ei23 Ld3)C3, (59) 
P"' = -Pa Ap-£i23L(dl Ad2)(Ci AC2), (60) 

^ = 7rp-pP,-a]J-ei23L(di Ad2)c3. (61) 

One can note that the dyadics e' and /Z"^ have a quite restricted form. In the case DC = 
they actually do not have an inverse, which is in contrast to the QDCM case. Actually, 
P-mcdia and Q-mcdia can be transformed to one another through Hodge duality [15]. 
The same property is also valid for the generalized Q- and P-media. 



3.2 SDCM 



Let us now consider the special case of the condition (44) simphfied by 7 = and (3 — 1, 
which yields the following first-order dyadic equation, 

aejv \}^^^^ + + MJ = AB + B A. (62) 

Equation (62) can be interpreted so that the symmetric part of the dyadic Mg — AB must 

be a multiple of eAr[I^^^^. Redefining a, the modified medium dyadic must thus be of the 

form _ _ 

Mg = aejv \}^^^^ + AB + A, (63) 

where A e E2E2 is an arbitrary antisymmetric dyadic. It is known that any antisymmetric 
dyadic mapping two-forms to bivectors can be expressed in terms of a trace-free dyadic 
Bo e EiFi as [21] 

t^er,lC\0of, tr%^0. (64) 
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It is now easy to verify that a medium defined by (63) satisfies the decomposition condition 
(42). In fact, any plane wave in such a medium satisfies 

# • * - a* • $ = $1 (e^ [M - ae^ [P)"^) I* 

= *|(AB + I)|* = (A|*)(B|*) = 0. (65) 

Alternatively, we can replace the solution (63) by 

Mg = aejv J^^^^ + I + AB + B A, (66) 

which corresponds to 

M= aP)^ + £jvL^ + eivL(AB + BA). (67) 

Without losing the generahty, the bi vectors can be assumed to satisfy A • B = 0, whence 
the last one of the three terms in (67) is trace free. In this case the three terms correspond 
to the respective axion, skewon and principal parts of the medium dyadic M. While the 
axion and skewon parts may be arbitrary, the principal part is restricted to be of the 
simple form as defined by the two bivectors A and B. Since the principal part is not 
complete, i.e., it does not have an inverse, some trouble in interpreting the medium in 
terms of three-dimensional medium dyadics may be expected. If A = B is chosen in 
(66) and (67), SDCM reduces to a simplified class of media, previously called that of 
doubly-skew media [22]. 

3.3 3D expansions for SDCM 

Because SDCM defines a novel class of decomposable media, it is interesting to find its 
definition in terms of three-dimensional medium parameters. Let us expand the trace-free 
dyadic Bp of (64) as 

Bo = Q-Fe475-FCs£4-e4£4(trCs). (68) 
where is a spatial dyadic, a spatial vector and 7^ a spatial one-form. Applying 

U T, + 6464, = f ^ + Ue4e4, (69) 

where is the spatial unit dyadic, we have 

= (Cs^hf - £4 A A C3 - 7s A A 64 - £4 A {t, - tits Tf ) A 64 (70) 
Further, we can expand 

A = ei23 [ocs + A 64, B = 6123 [f^s + A 64, (71) 

where the vectors a.s,hs and the one-forms cxs,/3s sue spatial. In terms of these we can 
write 

£jv L( AB + B A) = (-£4 A + £123 La3) (6123 [f3, + b, A 64) 
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+(-£4 A /3, + £i23Lbs)(ei23L«5 + a, A 64). (72) 

Inserting the expansions in (67) and equating with (21) we can identify one set of three- 
dimensional medium dyadics as 

^ = aTfT ^ (^^AQT + s,23ei23 LL(as/33 + b,a,), (73) 

e = -7. A ]J + £123 L(a.b, + b,a,) , (74) 
= -TJ Ac, - (a,/3, + /3,a,)Jei23, (75) 
I = -aTf - {ts - trZJsf - {(Xshs + (3,Sis)- (76) 

It can be seen that the dyadics 6123 [e' and JI^^\si23 have arbitrary antisymmetric 
parts while their symmetric parts are incomplete consisting of only two dyads. The 
decomposition (42) can now be written for the 3D fields as ^ 

(B|(ei23La.) +E|a,)(B|(ei23L/3J + E|b,) - 0. (77) 

The 3D Gibbsian-dyadic representation corresponding to (73) - (76) can be obtained 
through the transformation rules (25) - (27). It turns out that their analytic expressions 
become quite extensive and they are omitted here. 



3.4 Example of SDCM 

As an example, let us consider SDCM with vanishing magnetoelectric parameters: 

|^ = J^ = 0, ^ ^ = 0, ^ = 0. (78) 

This implies 

a = 0, ts = 0, a,/3, + b,a, = 0, (79) 

in the above 3D representations. Limiting to the most general case, i.e., that none of the 
quantities a,, /3^, b,, cxg vanishes, we must have b, = Aa, and /3g = — Aa, for some scalar 
A. Thus, in this case we can write 

e = -Is A t + 2A(£i23 La.)a„ (80) 

= -If A c, + 2Aa,(a,Jei23), (81) 
and the corresponding Gibbsian dyadics become 

^g= ^^^J^{CsCs + 2X{Cs\0Cs)ei23[i0CsA]'^)), (82) 

f5 = 2Aa,a,-ei23L(7.Af). (83) 

The medium defined by these expressions is characterized by both electric and magnetic 
gyrotropy. For example, the permittivity dyadic (83) can be approximately realized by 
magnetoplasma in a high static magnetic field for low frequencies [23] . 



^Here B stands for the magnetic two-form and not for the bivector. 
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4 Dispersion equations 

To verify the decomposition of plane waves let us derive the dispersion equation for the 
general plane wave in all of the previous medium cases. The equation for the potential 
one- form can be obtained by starting from (33), (34) which yield 

e^[(i/ A *) = -r/J(e^L*) = -(ryjM,H|(/> = 0. (84) 
This can be expressed as 

B{u)\(f> = 0, U(i/) = uu\\Mg, (85) 

where D(i/) e EiEi is the dispersion dyadic. The axion part of M^, does not contribute 
because 

i/i/JJ(ejvLT^^^^) =ejv[(i/Ai/ aF) = 0, (86) 
and we can omit it in all medium cases. One may note that (85) implies 

(u) [(f> = 1 0) A (u) = 0, (87) 

which will be applied in the sequel. 

4.1 QDCM 

In the case of DCM the dispersion dyadic equals that of the generalized Q-medium, 

B(i/) = ivi/J J (M^(2) + DC) . (88) 

This case has been analyzed in [19], but let us retrace the steps for convenience. For 
simplicity, let us define the vectors 

d = ivjD, c^u\C. (89) 

Expanding 

and applying the rules valid for normalized Q, 

= {lyiymW''^ = eive^vLL^^'^, (91) 
(i/i/JJ^(2))(3) ^ (^^11^)2^(4) ^ (i,i,||^)2e^e^r, (92) 
(90) can be expanded as 

+(e^^e^LL(^"'^>^^)Adc) 

= M\eNeN[[r^u){uiy\\Q){M{iyiy\{Q)+^-'^\\dc). (93) 

Equation (87) yields a scalar dispersion relation which splits in two quadratic equations 
as _ 

u\Q\u = 0, (94) 

i/|(M^ + CD[[^-^)|i/ = 0. (95) 

Actually, (94) corresponds to the A-wave and (95) to the B-wave as will be shown in the 
next section. 
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4.2 PDCM 

Neglecting again the axion term, PDCM coincides with the generahzed P-medium whose 
dispersion equation was derived in [15]. Omitting the details, quite similar to those of the 
QDCM, the dispersion equation can be split in two equations which are of the form 

u\{Blf)\u = Q, (96) 

u\{C[p-')\u = 0. (97) 

4.3 SDCM 

For SDCM the dispersion dyadic equals 

U(i/) = i/iyJJ(I + AB + BA). (98) 

Expressing the general antisymmetric dyadic as in (64), where Bo e EiFi may be any 
trace-free dyadic, we can write 

^^Jjl = -u\{e^l{%^]flu) = FLF, (99) 

where the bivector defined by 

F^BnUuAu'), u' = W,\u, (100) 

is simple since it satisfies F • F = 0. The dyadic F [F e EiEi is antisymmetric and it can 
be shown to satisfy 

(FLF) (2) ^ FF, (F[F)(^) = 0. (101) 
Defining for simplicity the vectors 

i/JA = a, i/jB = b, (102) 

we can expand 

W\v) = (F[F + ab + ba))(') 
= (FF)^(ab + ba) + (FLF);^(ab + ba)^^). (103) 
Applying the orthogonality i/|a = u\h = 0, we obtain 

F Aa= {eN[{t^ Au')) A3i^ -(ejvLi')(i^'|a)> (104) 
and similarly for F A b, whence 

(FF);^(ab + ba) = 2{eNeN[[j^J^){i^'\a){u'\h). (105) 

The last term of (103) vanishes due to 

(F[F);^(ab;^ba) = {eN[tyy\{ah - ba) A (a A b) = 0. (106) 
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Thus, we are left with 

W^\u) = 2{eNeNlluu){u'\a){u'\h), (107) 

which vanishes due to (87). Thus, we must have either i/'|a = or i/'|b = satisfied by 
u'. 

In conclusion, in the SDCM case, the fourth-order dispersion equation splits in two second- 
order equations 

u\{%\A)\u = 0, u\{%\B)\u = 0. (108) 

The plane-wave propagation depends on the metric dyadics BqJ A and BqJ B belonging to 
the space EiEi. The medium has no birefringence if the symmetric parts of these two 
dyadics are multiples of one another. When in addition A is a multiple of B, the medium 
coincides with the doubly- skew medium of [22]. 

5 Properties of the plane-wave fields 

Let us now check whether the plane-wave fields satisfy the decomposition properties as- 
sociated to the corresponding media. 

5.1 QDCM 

Applying (50) we can write for the QDCM case 

B{i^)\(f> = i/ivJJ(M^(2) + DC - /5eArLl^')^)|0 

= M(i/j/||Q)Q|0- M(Q|i/)i/|Q|</) + dc|0 = 0. (109) 

Since the potential one-form is not unique, one can choose an additional condition (Lorenz 
condition) without changing the field two-form $ = i/ A 0. Choosing the condition 

i/|^|</) = 0, (110) 

(109) is reduced to 

M(i/zy||^)^|</) + dc|(/) = 0. (Ill) 
Now the solution corresponding to the dispersion condition (94) leads to 

c|</>= -C|$ = -A|$ = 0, (112) 

because of the definition (47). Thus, (94) corresponds to the A- wave. 

To verify that the field corresponding to the dispersion equation (95) satisfies the B-wave 

condition B|# = appears more complicated. Starting from c\<p ^ and i^i^HQ 7^ 
while assuming the condition (110), from (111) we see that the potential must be of the 
form _ 

= AQ-^|d, (113) 
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where A is some scalar coefficient. The condition (110) then requires that v satisfy 
f|d = 0. Invoking (52) and the rule 

2D A(i/JD) = i/J(DAD) = (ejvM(D-D), (114) 

after a few algebraic steps we obtain 

B|* = AB|(i/A^-^|d) 

= -(D-D)(M(iy|^|i/) + c|^-^|d). (115) 
Comparing with (95), the corresponding solution can be identified as the B-wave. 

5.2 PDCM 

Omitting again some details, in [15] it has been shown that the field corresponding to the 
dispersion equation (97) satisfies the condition 

C|$ = 0, (116) 

i.e., it represents an A-wave. Similarly, the solutions of the dispersion equation (96) 
correspond to the field condition 

(MDp2) + l(D.D)C)|$ = 0, (117) 

which equals B|$ = due to the definition (55). 

5.3 SDCM 

For the SDCM case, the equation for the potential one-form (85) with expressions from 
Section JVC inserted and the axion term omitted, yields 

n(i/)|0 = (F[F + ab + ba)|(/) 

= ejvL(i^ A I/' A 0) + a(b|</)) + b(a|(/)) = 0. (118) 
Multiplying by 0| yields (a|(/>)(b|0) = 0. Multiplying by v'\ leads to the condition 

(i/»(b|<^) + (r.'|b)(a|<^)=0, (119) 

from which we obtain the relations 

i/'|a = i/|(IoJA)|i/ = ^ a|0 = -A|* = O, (120) 

iv'|b = i/|(l<,jB)|iy = ^ b|0 = -B|$ = O. (121) 

This shows us that the plane waves corresponding to the two dispersion equations (108) 
are respectively A- waves and B- waves. 
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6 Conclusion 



In this paper we have apphed four-dimensional differential-form formalism to define media 
in which field two-forms can be decomposed in two parts: the A-fields and the B-fields: 
$ = + ^B- The decomposed fields are defined in terms of two given bivectors A 
and B so that they satisfy $a|A = and $b|B = 0. Media with such a property 
arc called decomposable media (DCM) for brevity. It is shown that a large class of 
four-dimensional DCM medium dyadics must satisfy an equation which is either linear or 
quadratic. The media satisfying the linear equation are called special decomposable media 
(SDCM). Those satisfying the quadratic equation are shown to fall in two subclasses which 
equal those previously known as generalized Q-media and generalized P-media with an 
added axion term. Such subclasses have been dubbed QDCM and PDCM, repectively. 
Basic properties of all three classes of decomposable media arc discussed in the paper. 
The properties of QDCM and PDCM are based on earlier studies on generalized Q- 
and P-media while the SDCM class appears to be a novel generalization of the class of 
non-birefringent doubly-skew media. Dispersion equations for a plane wave are derived 
for all classes of media and their correspondence to the decomposition of the field are 
demonstrated. In the Appendix it is briefly shown that a medium dyadic satisfying the 
quadratic equation • M = cuejv [T^^^^ must belong to either the class of Q-media or that 
of P-media. 



Appendix: Quadratic dyadic equation 

Let us consider the following quadratic equation for the medium dyadic M e F2E2, 

. M = M^|(ejv[M) = ae^ (122) 

in the case a 7^ 0. Inserting the three-dimensional expansion (21), (122) equals the set of 
four 3D conditions 

^"'^|(ei23L^) +^|(ei23L^"') = 0, (123) 
P~'^|(ei23L^) +^|(ei23L^) = -«ei23Lf , (124) 
F|(ei23L^) +f'^|(ei23L^-') = -«ljei23, (125) 
F|(ei23Le)+f^1(ei23L^) = 0. (126) 
Since (124) equals (125) transposed, we can ignore (125). 

Now one can show that, out of the two pairs of dyadics a, JT^ and e', /3, exactly one dyadic 
in each pair possesses a 3D inverse. To prove this, let us first consider the pair a,]I~^ and 
assume that neither of the dyadics possesses an inverse. Choosing the reciprocal bases 
ej, £j in a suitable manner, we can expand 

T^}eu3 = SiKi + £2^2 e FiFi, (127) 

ei23ei23LL^ = eiai-Fe2a2 + e3a3 e EiFi, (128) 
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where Ki, K2 and ai,a2, 0:3 are one- forms of which the latter three are hnearly dependent. 
These substituted in (123) yields a dyadic equation between the one- forms 

Kiai + K20t2 ^ -{0CiKi + a2K2). (129) 

Assuming that cki A 0:2 7^ 0, (129) implies that ki and occupy the same subspace as 
oci and CX2, and that furthermore, for some scalar A, we have 

Ki — Act2, K2 — —Acti. (130) 

The relation now becomes 

^"^ = A{£iOi2 - e2Q:i)Jei23 = -AeaJ^. (131) 
Substituting this in (124), the resulting equation 

^\{A(\f)[e,)\{e,23\j)+e,23m = -«ei23L? (132) 

does not have any solutions since the left-hand side has no inverse while the right-hand 
side has one for a 7^ 0. Thus, the assumption that neither a nor fl'^ have an inverse is 
incorrect. Releasing the constraint ai Aa2 7^ 0, the case ]I~^\£i23 = si^, ei23ei23LL^ — 
can be shown to lead to the same result. Thus, at least one dyadic of the pair a, 'pT^ must 
possess an inverse. The same conclusion is valid for the pair f', (3. 

To demonstrate that exactly one dyadic of each pair possesses an inverse, we notice that 
the left-hand sides of (123) and (126) are symmetric dyadics, whence all four terms must 

be antisymmetric dyadics. Because 3D antisymmetric dyadics A can be expressed in terms 
of some vector a or one-form a as 

leEiEi, ^ I = ei23L(aATf), (133) 

AGE2E2, ^ I = aALjei23, (134) 
the dyadics in (123) and (126) must have the general form 

^|(ei23L^~') =aAT,Jei23, (135) 

F|(ei23L^')=ei23L(«ATJ). (136) 

Since antisymmetric 3D dyadics do not have an inverse, neither do the left-hand sides 
of (135) and (136). This concludes the proof. As a conclusion, there are exactly two 
invertible dyadics among the four dyadics. Let us split the problem of finding solutions 
to (122) by considering the four possible cases separately. 

6.1 Case 1: (a, 13) 

Assuming that a and (3 possess inverses, the other two dyadics can be solved from (135) 
and (136) as 

= (£1236123 LL^~'^) A a, (137) 
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e= (£1236123 LL^~'^)L«- (138) 

Defining the dyadic 

^= (£1236123 LL^)I^ (139) 

the inverse of which exists due to the above assumption, we can express (124) as a con- 
dition for the dyadic X, 

£1236123 L [(^"'"^Aaa) + X= -alj. (140) 
Now it is easy to check that the solution of (140) must be of the uniaxial form 

X=ATf + ^aa, detf^A{A^ + Si\a), (141) 

where the parameter A and the scalar a of (122) are related as 

a = - ^ (^^ + a| a) = ^detX. ( 142) 

At this point wc can construct the medium dyadic M satisfying (122) for Case 1 from 
(139), (137) and (138) as 

M = ?-^(detp(2) _ A^^0\cx + det^£4)(a + ^64)), (143) 

^Met^ 

and it can be further expressed in the compact form 

M = mP^), M= L — =, (144) 

AMetXdet/^ 

P = (detX)^ - A0\a + det^£4)(a + ^64). (145) 

Media defined by medium dyadics of the form MP*^^) for some dyadic P have been called 
P-media [15]. 

6.2 Case 2: (1',^"^) 

Assuming that t and Jl~^ possess inverses, the other two medium dyadics can be expressed 
as 

^ = -(£1236123 LL^^) A a, (146) 

^=-(£l236l23LL^~''^)L«- (147) 

Defining the dyadic 

^= (£1236123 LL^"'^) I e eFiEi, (148) 

(124) can be shown to yield (140) in spite of the different definition of X. Thus, the 
expression (141) can be applied. Let us for simplicity define the additional dyadic 

= _ =(3) 

Y = ei23L7^ eEiEi, Ay = £i23£l23||Y (149) 



18 



in terms of which we can express from (148), (146) and (147) the dyadics 

t = eu3[^^\% (150) 

^=-ei23L^^Aa, (151) 

^= ^ei23L^(^)nX(2H«- (152) 

AydetX 

Substituting these in (21), the modified medium dyadic M^, = eAr[M becomes 

% = + - ^4)(^ - 64). (153) 

This can be expressed in the compact form 

^ = 7^|X- ^{W\cx- A^AYe^){a- Ae^). (155) 

Media defined by modified medium dyadics of the form = MQ^^^ for some dyadic Q 
have been called Q- media [14] . 

6.3 Case 3: {a,f) 

As a third case we consider the possibility that a and e have inverses. Equation (124) 
now becomes 

a A (ei23L^"^|f') + (^l^'"^^Jei23)L« = "6123 Jf- (156) 

Denoting this time 

X = ei23L^-'|e eEiEi, = £iveiv||^^'\ (157) 

which has an inverse within the limits of Case 3, the condition (156) takes the form 

a A X + (61236123 LL^~'^)L« = «6i23Lt- (158) 
Multiplying by |a we obtain 

a A X|a = Q;ei23[Q:, a|a = 0, (159) 

whence (158) becomes 

bAX=aei23jf, b = a+^. (160) 

Because of the assumption a 7^ 0, this condition leads to an impasse: the dyadic on the 
right-hand side has an inverse while that on the left-hand side has not. Thus, the original 
assumption that a and f' possess inverses is obviously invalid. 
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6.4 Case 4: i(3,JI-^) 

As the final case we consider the possibihty that j3 and /x"^ have inverses. Equation (124) 
now becomes _ _ _ 

(^-i^|r'^Jei23)L« + aA (ei23L^i;5) = -ae,23\J^ . (161) 

Since this is of the same form as (156), the conclusion must be similar: the assumption 
that /3 and have inverses cannot be valid for M to satisfy (122) with a 7^ 0. 

In conclusion, there are only two classes of solutions to the quadratic equation (122), that 
of Q-media and that of P-media. 
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